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We construct and discuss a toy model of the population of numerous non-identical extragalactic 
sources of ultra-high-energy cosmic rays. In the model, cosmic-ray particles are accelerated in 
magnetospheres of supermassive black holes in galactic nuclei, the key parameter of acceleration 
being the black-hole mass. We use astrophysical data on the redshift-dependent black-hole mass 
function to describe the population of these cosmic-ray accelerators, from weak to powerful, and 
confront the model with cosmic-ray data. 
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I. INTRODUCTION 

The origin of ultra-high-energy (UHE) cosmic rays 
(CRs), that is of cosmic particles with energies E > 
10 19 eV, is presently unknown. However, there are nu- 
merous hints, both in data and in theory, which might 
help to constrain possible models of UHECR sources. In 
particular, the observation [lH3{ of the spectrum steepen- 
ing consistent with the Greisen-Zatsepin-Kuzmin [j, Q 
(GZK) cutoff, together with the global isotropy of the 
arrival directions (see e.g. @) and the fact that the 
UHE particles are not expected to be confined by the 
Milky- Way magnetic field (see e.g. [7|) suggest that the 
bulk of cosmic rays at these energies have extragalac- 
tic origirQ. Next, the lack of clustering of arrival direc- 
tions at small scales is a powerful tool [9( to constrain 
the number density of sources. In recent data of the 
Pierre Auger Observatory (PAO), no evidence for clus- 
tering at E > 5 x 10 19 eV is seen [l(| which translates 
into the number density of n > 1CU 4 sources per cubic 
Megaparsec, which means that sources of these extreme 
particles should not be exceptional and, most probably, 
some of them should be located relatively nearby. The 
latter fact gets further, though limited by statistical sig- 
nificance, support from the shape of the GZK feature in 
the spectrum which does not seem to be very sharp (cf. 
Ref. [13). 

On the other hand, it is a nontrivial task to find partic- 
ular astrophysical objects which could serve as UHECR 
accelerators. Even without a detailed modelling of the 
acceleration process, a number of simple estimates rule 
out many classes of potential sources. These simple crite- 
ria include in particular the geometrical (Hillas) criterion 
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fl2j and estimates of radiative energy losses of particles 
being accelerated (see e.g. Refs. [lj, EH)- Analysis of 
the modern astrophysical data demonstrates [l5[ that the 
combination of these constraints leaves just a few candi- 
date classes of sources capable of acceleration of particles 
to UHE energies. Leaving aside large-scale structures 
where interaction losses are expected to suppress the en- 
ergy gain, the conventional diffusive (e.g., relativistic or 
non-relativistic shock) acceleration may work only in ul- 
trarelativistic jets, hot spots and lobes of exceptional ac- 
tive galaxies (powerful radio galaxies and blazars) which 
are not that abundant in the nearby Universe. For very 
special field configurations when synchrotron losses are 
suppressed and the curvature radiation dominates, pos- 
sible acceleration sites include also gamma-ray bursts 
(GRBs) and immediate neighbourhood of supermassive 
black holes (SMBHs) in the galactic nuclei. While it is 
unclear whether these field configurations may be present 
in GRBs, recent IceCube results disfavour the GRB sce- 
nario anyway [l6[ (see however Ref. [13) ■ At this level 
of reasoning, the SMBH environment remains a viable 
option. 

A natural assumption is that numerous UHECR 
sources are not identical - there should be less and more 
powerful accelerators where the maximal energies, in- 
jection spectra and fluxes of accelerated particles are 
different. Until now, numerous attempts to model the 
sources of UHECRs and to confront theoretical predic- 
tions with experiments often assumed that these param- 
eters are fixed once and for the entire Universe (see e.g. 
Refs. fl8l . [Tij and numerous other works; see however 
Ref. [2Cf where acceleration in non-identical jets was con- 
sidered). While, for numerous sources, the assumption of 
equal fluxes is well justifiable (in the sence that only the 
mean flux of a large sample of sources is important and 
this mean flux does not vary significantly from one region 
in the Universe to another) and the injection spectrum is 
often fixed by the acceleration model, the maximal ener- 
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gies are expected to vary significantly. As it was recog- 
nised in Ref. [2l|, these variations affect the observable 
spectrum seriously. In this work, we attempt to present a 
toy model of numerous and different sources of UHECRs 
and, within certain assumptions, to confront it with the 
experimental data. 

To this end, we choose a simple toy model of particle 
acceleration in the immediate vicinity of SMBH put for- 
ward in Refs. [22|, (HI- The reason to choose this particu- 
lar model is twofold. First, unlike many other models, it 
allows 15] for UHECR acceleration in numerous nearby 
sources. Second, as we will see below, within some realis- 
tic assumptions, the acceleration capabilities of a source 
are determined by a single parameter, the SMBH mass. 
At the same time, the demography of SMBHs is well 
studied by astrophysicists and we take this advantage to 
describe the population of sources easily. 

The rest of the paper is organized as follows. In 
Sec. HH we give a brief review of the acceleration model 
of Refs. (22I |23|| and make a bridge between the parame- 
ters which determine the maximal energy of accelerated 
particles and the SMBH mass. In Sec JIIH we discuss 
the astrophysical data on the SMBH population, merge 
them with the acceleration model, calculate the spectrum 
of UHECRs with the account of propagation from source 
to the observer and compare it with the experimental 
cosmic-ray data. We obtain a good agreement with the 
observed spectrum by fitting the spectrum with only two 
continuous parameters, the overall normalization and a 
single free parameter of the model. Sec. IIVI demonstrates 
that the population model of Sec. IIIII satisfies simple ob- 
servational constraints: it does not produce too much 
secondary gamma rays, it results in an acceptable num- 
ber density of sources and, with the best-fit normaliza- 
tion, it does not require enormous luminosity of a single 
source. We give our conclusions and discuss our results 
in Sec. M 



near the rotation axis, particles moving along magnetic 
lines are accelerated by the electric field. 

In this case, radial components of the electric and mag- 
netic fields in units h = c = G = 1 in locally non-rotating 
frame in Boyer-Lindquist coordinates on the symmetry 
axis are: 
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where M is the black hole's mass, a < M is its angular 
momentum per unit mass, r is the radial coordinate and 
Bq is the external homogeneous magnetic field. 

Neglecting for the moment the energy losses, the maxi- 
mal energy gain of the accelerating particle with a charge 
Ze is determined by the available potential difference in 
the region along the rotation axis, 
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Thai — M + %/ M 2 — a 2 is the radius of the black-hole 
horizon where particle acceleration starts while r max lim- 
its the size of the region along the rotation axis, where 
acceleration is possible (in Ref. [23| , it is called "the vac- 
uum gap" due to the absence of numerous charged parti- 
cles in this region except of the single test particle being 
accelerated, which does not change the electromagnetic 
field). 

Of course, one can rewrite Eq. ([2]) for the potential dif- 
ference and £ max in terms of the distance-averaged elec- 
tric field E, 

£ max (a) = ZeE(a)H, 



II. A TOY MODEL OF PARTICLE 
ACCELERATION IN THE BLACK-HOLE 
MAGNETOSPHERE 

A toy model of particle acceleration in the black-hole 
magnetosphere was proposed by Neronov et al. [22I [23|. 
Let us briefly discuss the model and its parameters. 

Assume that a stationary rotating black hole without 
electric charge is embedded into the external magnetic 
field, homogeneous at the horizon distance scale. In gen- 
eral, the magnetic field is inclined at some angle \ with 
respect to the black-hole rotation axis. There is a well- 
known exact solution of Maxwell's equations in the Kerr 
metric for each inclination an gle y of an asymptotically 
homogeneous magnetic field [2J, [25|. For instance, if 
X = 0, then a rotation- induced electric field is parallel 
to the magnetic one on the symmetry axis and its direc- 
tion depends on the directions of both the magnetic field 
and the black hole's rotation velocity. Thus in the region 
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where H = r hor - r max . 

By making use of Eq. (pj one obtains 
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It is easy to see that the expression in square brackets 
equals to —H(H + 2M r^or — 2a 2 ) and therefore is always 
negative for any value of H (remember that the angular 
momentum per unit mass a varies between zero and M ). 
Thus for parallel magnetic field and angular momentum, 
aM > and the radial component of the electric field Ep 
on the black-hole rotation axis is negative. So if there are 
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maximal energy of an accelerated particle is 



FIG. 1. Dependence of the normalized maximal energy £ from 
the angular momentum a of the black hole with mass M for 
different values of the vacuum-gap size H (curves labeled by 
values of H/M; see the text for definitions of parameters). 



negative charges near the rotation axis, they will be ac- 
celerated away from the black hole. According to Eq. 
the difference between their energies at r max and r/,. or 
is positive, so they gain energy while moving along the 
rotation axis away from the black hole. In the opposite 
case, when the magnetic field and angular momentum are 
antiparallel, the radial component of the electric field on 
the black hole rotation axis is positive. So positive-charge 
particles, situated near the rotation axis, are accelerated 
away from the black hole; according to Eq. ([3]), the en- 
ergy gain in this case is positive for the positive-charge 
particles. 

In Ref. 23], a simple expression £ max ~ ZeB H 
rather than Eq. ([3]) was used; however, more precisely, 
Eq. (J3J implies £ max < ZeB^H. The dependence of 
£ = fmax/ ZeBoH from the angular momentum a for dif- 
ferent H is shown in Fig. [TJ For black holes with angular 
momentum a > 0.1M and H ~ (1 - 2)R S = (2 - 4)M, 
we have £ max ~ 0.1 ZeB H; however, for slowly rotat- 
ing black holes, a < 0.1M, £ max varies from zero to 
~ O.OlZeBoH . Hereafter, R§ denotes the Schwarzschild 
radius. 

However, the precise value of £ max is often irrelevant 
since particles cannot achieve this maximal energy be- 
cause of inevitable energy losses associated with the ac- 
celerated motion of the particle. The particle energy is 
determined by the balance between the energy losses and 
the energy gain per unit time, 
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It was shown in Ref. [23[ that protons can be accelerated 
to the energies of about 10 20 eV only if the magnetic 
field is almost aligned with the rotation axis. In this case, 
only the curvature radiation is relevant for an accelerated 
particle, and not the synchrotron one. We should also 
note here that we do not consider energy losses related 
to interactions of accelerated particles in the source. The 
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where B is the external magnetic field, R ~ Rs/x ^ s 
the curvature radius of magnetic-field lines, Ze is the 
particle charge, Am is the particle mass (A is the atomic 
number and m is the nucleon mass) and n is a coefficient 
between the electric field and the external magnetic field, 
\Ep\ = kBq (see Eq. (JTJ). Note that k is a function of the 
angular momentum a and the coordinate r. We point 
out that numerically, £ cutv < f max , cf. Ref. [l5j], and 
therefore it is £ curv which determines the actual maximal 
energy. For simplicity, we assume that all of the particles 
start with equal initial conditions and so all of them are 
accelerated to the same energy £ curv . 

The original model of Ref. [23| treated the magnetic 
field Bq as a free parameter. However, one may note 
that the field is constrained and, in particular, can- 
not be too high (see Ref. 
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sion). The maximal value of the magnetic field is deter- 
mined [26|, [2?J by the so-called Eddington limit, B^d = 

/ \ — 1/2 

( wHf . ) ^ Though this estimate may be ob- 
tained in several different ways, its most transparent 
meaning is that the magnetic-field energy density is equal 
to that of the accreting plasma, corresponding to the 
Eddington luminosity. To obtain the maximal possible 
particle energy, we should assume that the external mag- 
netic field is equal to the Eddington limit, Bq = i?Ed- 
However, below, we will find the spectrum of cosmic rays 
in the frameworks of this model. For this purpose we 
are interested in the actual maximal particle energy as a 
function of black-hole mass M rather than in the upper 
limit. So we should recognize that the realistic magnetic 
field can differ from the Eddington limit. In a general 
case, we can parametrize the external magnetic field as 
follows: 



B = kB 
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where a and k are some parameters. We note in pass- 
ing that several realistic models predict this kind of 
dependence, e.g. the Shakura-Syunyaev model [H, [29| 
(k « 0.31, a = 0) or the model of Ref. [30] (Jfe w 0.0093, 
a m —0.31; see Fig. 2 of Ref. [l5[ for comparison with 
scarce observational data). We will consider these op- 
tions below. We have 
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During acceleration, the particle emits curvature pho- 
tons. In what follows, we will need to obtain an upper 
bound on this emission. The peak energy of the photons 
is determined by the particle energy £ , 

f - 3 g3 
7 2m 3 R' 

and for the upper limit, we take £ — £ C urv, 
A 3 ( M 



~ 14 TeV , . 
7 Z 3 / 4 Vl0 9 M o 



The ratio of luminosities in photons, L 7 , and in cos- 
mic rays, Lqr, may be estimated by comparing the total 
available potential difference in the acceleration region 
along the rotation axis to its fraction, spent on the par- 
ticle acceleration: 
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where (C^ 1/4 ) ~ (0.1 - 2) for < a < M and R s < 
H < 6i?s- In numerical calculations presented below, we 
use k 1 / 4 = 0.7 and £ = 0.25, cf. Figs. EJ [3] 

To summarize, the model we use assumes a monochro- 
matic spectrum of accelerated particles with £ — £ C urv, 
Eq. (HJ) , in each particular source. The value of £ curv de- 
pends, within the assumed magnetic-held model, on the 
SMBH mass M only (in what follows, we do not consider 
acceleration of other particles than protons). The overall 
flux from the source remains a free parameter. 



III. POPULATION OF THE SOURCES AND 
THE OBSERVED SPECTRUM 

Having discussed the model of particle acceleration in 
a particular source, we switch now to the population 
of sources. As we have seen, the properties of a single 
source are determined, within the magnetic-field model 
we choose to study, by the SMBH mass (the dependence 
from the spin is weak). To reconstruct the UHECR spec- 
trum one has to consider the population of SMBHs dis- 
tributed in mass and luminosity. For simplicity we will 
assume that the mean SMBH luminosity in cosmic rays 
is related to its mass, 



L C r oc M 13 , 
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where /3 is an additional model parameter. Note that 
not every black hole can work as a source, because the 
source should possess some special properties (for exam- 
ple, small inclination angle, the vacuum gap larger than 
i?s, absence of numerous charged particles in the vicin- 
ity of the black hole which might imply a thin or even 



absent accretion disc). The fraction of sources where the 
mechanism works is also encoded in the mean luminos- 
ity, Eq. ([7]). The observed spectrum can be obtained by 
convolving the SMBH mass function with the (monochro- 
matic) single-source spectrum and Eq. (|7j and taking into 
account the propagation effects. 



A. The SMBH mass function 

Since dynamical measurements of SMBH masses are 
available for a very limited number of cases only, it is 
a difficult task to find the mass function precisely. For- 
tunately, a number of indirect methods to estimate the 
SMBH mass are available (see e.g. Ref. [3l| for a review). 
Despite having large uncertainties in individual measure- 
ments, these methods are suitable for obtaining average 
characteristics of the SMBH population, see e.g. Ref. 32 1 
for a review. 

An important feature of the SMBH mass function is 
its evolution: SMBHs grow fast, cf. e.g. Ref. 33]. For 
our calculation, we use one of the most recent published 
rcdshift-dependent mass functions [lij . Of two functions 
presented there, we choose to use the one based on the 
stellar mass functions because it has smaller statistical 
uncertainties. The systematic uncertainties of the mass 
function may be judged from Ref. [32[ and are well within 
the overall precision of our toy model. 



B. The observed spectrum 

Before reaching the Earth, the accelerated protons may 
interact with the cosmic microwave background. The 
main two processes modifying the shape of the propa- 
gated cosmic-ray spectrum are photopion production and 
e + e~-pair production. The former leads to a strong sup- 
pression of the proton flux above few tens of EeV known 
as the GZK effect [J, why the latter mostly domi- 
nates the attenuation below 10 EeV leading to the so- 
called "dip" feature in the spectrum l35l . |36| . We use the 
numerical code developed in Refs. [33]. The code also 
traces secondary particles produced in the interactions. 
It makes use of the kinetic-equation approach and cal- 
culates the propagation of nucleons, stable leptons and 
photons using the standard dominant processes (see e.g. 
Ref. [H). 

Fig. [U presents the predicted cosmic-ray fluxes in the 
best-fit model for the Auger spectrum [39j for different 
dependencies of the magnetic field Bq on SMBH mass 
mentioned in previous section. The red curves corre- 
spond to Bq given by the Eddington limit, t he g reen 
one describes the Shakura-Syunyaev model 1281. 1 2911 and 
the blue one corresponds to the model of Ref. [30] . The 
overall flux normalization is a free parameter. Besides 
we tried two values of angle x = 1° an d 5° and varied 
the luminosity dependence ([7]) parameter f3 in the range 
— 1 < f3 < 2. The best- fit parameter values are indi- 
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FIG. 2. The cosmic-ray flux predicted by the model for three 
different assumptions about the SMBH magnetic field, see 
labels on the plot and explanations in the text, versus the 
Auger experimental data [39j |. 



cated on the plot. One can see that the first two models 
produce satisfactory spectral fits above 10 EeV. 



IV. CONSTRAINTS 

In this section, we discuss additional consistency 
checks of the model. They include estimates of the 
accompanying gamma radiation which should not be 
in conflict with the measured diffuse gamma-ray back- 
ground, estimates of the concentration of sources and of 
the luminosity of a single source. We will see that the 
model passes these tests. For order-of-magnitude esti- 
mates in this section, we assume Bq ~ i?Ed- 



A. Concentration of the sources 

Let us check that the local concentration of sources 
of cosmic rays with energies f> 6 X 10 19 eV is not in 
conflict with the lower limit [Hj based on the statistics 
of clustering. To this end, we integrate the SMBH mass 
function over the range of masses corresponding to these 
energies. 

The dependence of the particle energy from the black- 
hole mass is given by Eq. (|5|). Every black- hole mass 
corresponds to a range of particle energies due to the 
variations in the value of k. The black-hole mass function 
includes all the black holes with fixed masses and so all 
the black holes with every value of n. Thus, the lower 
limit of the required mass interval is determined by £ = 
6 x 10 19 eV and the maximum value of k. It was shown in 
Ref. [23| that particles can be accelerated to the energies 
of about 10 20 eV only if the size of the vacuum gap is 
not smaller than the Schwarzschild radius. In this case, 
the highest value of k 1 / 4 is ~ 0.7 (see Figl3]). Black-hole 
masses corresponding to particle energies £ > 6 x 10 19 eV 
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FIG. 3. Dependence of the coefficient k 1 / 4 from the angular 
momentum a of the black hole with mass M for different 
values of the vacuum-gap size H (curves labeled by values of 
H/M; see the text for definitions of parameters). 



are then M > 10 7 Mq. Integrating the mass function in 
this range of masses, we obtain 



dn 



d log M 



d log M. 



The integral is saturated at its lower limit that is M max 
can be taken arbitrary high to obtain the following esti- 
mate, 
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As we see, the total concentration of sources is larger 



than the clustering lower bound of 10 



Mpc a 



As we 



have discussed above, only a fraction of the calculated 
concentration n corresponds to the true concentration 
of the sources. Our estimate tells us that this fraction 
should be not less than a few per cent which is reasonable. 



B. Luminosity of a single source 

A simple estimate of the luminosity of a single source 
may be obtained as follows. Consider the observed flux 
of cosmic rays with energies £ > 6 x 10 19 eV. The value 
j(£ ), which is often reported, is 
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where N(£) is the number of particles with energy < £ , 
per unit area per unit time. Then, for the flux F we have 



1 dF 
4^~d£' 



and, from the recent data [1, H| , we obtain the estimate 
F ~ 1.3 x 10 5 
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On the other hand, this flux is produced by the sources 
situated in the GZK sphere, where the GZK horizon ra- 
dius is i?GZK ~ 130 Mpc for energies £ > 6 x 10 19 eV, 
see e.g. [4(|. We assume that all these sources have ap- 
proximately equal cosmic-ray luminosities Lq \^f] , inde- 
pendent from the source black-hole mass (in reasonable 
agreement with the best-fit values of /?, Sec. IIIip . Be- 
cause in this case the GZK radius corresponds to very 
small redshifts z < 0.1, we can neglect changing of the 
source concentration with the redshift. For the order- 
of-magnitude estimate, we neglect also the difference be- 
tween the energy with which the particle was emitted 
and the final particle's energy with which we detect it 
on the Earth. Let us also note that, because we are in- 
terested in all sources of the accelerated particles with 
energies £ > 6 x 10 19 eV, in the expression for the to- 
tal flux we have to substitute the total concentration of 
sources, Eq. ©, that was obtained by integrating the 
mass function. 

The flux from every single source, situated at the dis- 
tance d from us, is Fq — Lq/ (4nd 2 ). Thus for the total 
flux we have 
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Using n ~ 10~ 3 T7 ^, we obtain L„~ 6 x 10 39 SE. The 
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corresponding luminosity in photons is L 7 = rj- Lq, where 



r\ is given by Eq. (J6)). One has 
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This value is much smaller than the typical bolomet- 
ric luminosity of an AGN, L AGN - (10 41 - 10 43 ) 
and by far does not exceed the Eddington limit, Led = 
10 47 ( 10 ^m o ) Taking the concentration of the 

■, does 
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sources of order of the lower limit, n ~ 10 
not result in a conflict as well. 

In all the cases, the luminosity is not that far from 
the luminosity of an AGN, and a natural question arises: 
can we see our sources in TeV as point sources? A simple 
estimate of the flux of TeV photons associated with the 
particle acceleration from one of the nearest sources, R ~ 
10 Mpc, gives the answer: we cannot see them as point 
sources because the flux from a single source is smaller 
than the sensitivity of the telescopes. Indeed, the flux 



of F 1 < 

io- 13 - 



10 



-1 eV 



corresponds to the counting rate of ~ 
beyond the reach of current TeV telescopes. 
Of course, this does not mean that strong TeV sources 
cannot accelerate UHE particles by this mechanism: the 
TeV emission may have a totally different origin. 



C. Diffuse gamma-ray background 

While individual cosmic-ray sources have quite low 
jamma-ray luminosities, one may wonder about the total 



emission of all sources in the Universe (beyond the GZK 
sphere). The emitted curvature photons have energies 
of order a few TeV and interact with the infrared back- 
ground radiation to produce electromagnetic cascades in 
which the energy of the leading gamma rays downgrade 
to the GeV band. Electrons in the cascade are deflected 
by cosmic magnetic fields so distant sources contribute 
to the diffuse gamma-ray background. Let us check that 
this contribution does not exceed the measured value of 
the diffuse flux. 

A simple estimate may be obtained as follows. Con- 
sider particles with energies £ ~ 10 20 eV. Cosmic rays 
with these energies arrive to us from the interior of the 
GZK sphere only, but the associated photons come from 
all sources at all distances. Knowing how much cosmic 
rays come from the interior of the GZK sphere, we can 
estimate how much of them are present in the Universe 
(keeping in mind that the number of sources depends on 
the distance). And supposing that the luminosity of a 
source in cosmic rays is connected with its luminosity in 
associated photons, Eq. ©, we can estimate the total 
photon emission from all the sources. 

The flux at the Earth, Eq. © , is expressed in terms of 
the energy at detection, £, which in general differs from 
the energy at injection, E m . In particular, the account 
of the Universe expansion (even neglecting additional en- 
ergy losses) results in £ = E in (l + z) -1 , where z is the 
redshift of the source. The number of the emitted par- 
ticles per unit time is also (1 + z) times higher than the 
number of detected particles on the Earth per unit time. 

The contribution to j(E) from the sources located at 
redshift z is 

1 rfn BH (M(£ in )) dN 1 

dj{E > z) = ~, dE~ ~dT W) 

where nBn(M(Ei n )) is the number density of black holes 
with masses < M(E- m ) (for M(E), see Eq. ©) at the 
redshift z, dNo/dt is a number of detected particles per 
unit time from one of the sources, which were emitted 
with the energy E- lni dV is the volume of a spherical 
layer at the distance z from us, S is the area of a sphere 
with the radius equal to the distance from the source 
to the Earth (it is necessary for calculating a number of 
particles through the unit area on the Earth). Taking 
into account the dependence from the redshift z, we may 
find the source cosmic-ray luminosity L{E- U1 ) 1 that is its 
total energy emission in cosmic rays per unit time, 

r (F ,dE t o t _ dN™ 

where dN^/dt is a number of particles with the energy 
E{ n emitted by the source per unit time and we assumed 
that the source emits particles with the only energy, de- 
termined by its mass M. 

For the flat Universe, we have S — 47ra§p 2 (z) and dV = 
Saodp(z), where 



Pi?) 



H(z)a 
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is the geodesic coordinate distance from the observer to the source location. Finally 
I 



m = 



1 dn BH (M((l + z)E))L((l + z)E) dlo &Tq 



G 



4tt 



dlog^ {1 + zfE 



dE w 



o 



a dp(z), 



(11) 



M((l+z)E)) 



where dnBH (M) / d log M is just the mass function for 
a given z [34j . z max is the redshift of the most distant 
source. Here we are interested in the values of the mass 
function at the points M((l + z)E) as a function of z. 
This function could be easily constructed using the data 
from Ref . [34] . Using Eq. ([5]) , we have 



o 



dE u 



(M((l + z)E)) 



1.16 

Ei n 



1.16 



(l + z)E 



For simplicity, let us suppose that all the sources have 
the same luminosities L((l + z)E) = Lq. This assump- 
tion is in reasonable agreement with the results obtained 
in section IIIII Thus we can take the luminosity out of 
the integral. Let us now imagine just for a moment, that 
cosmic rays with such energies could come from all dis- 
tances, and calculate the ratio of the values j{E) for the 
cosmic rays with the energies ~ 10 20 eV from the GZK 
sphere (z max = ^gzk ~ 0.01) and from the whole Uni- 
verse (e.g. z max = 2). The luminosities in front of these 
two integrals are cancelled and after integrating we ob- 
tain: 



E 2 ]{E) tot 
E 2 j(E) GZK 



20. 



Taking into account the observed value of the cosmic-ray 
flux, we estimate 



E 2 j{Ef 



10 5 



eV 



s • sr 



For calculating the value of E 2 j(E) 7 for the gamma ra- 
diation from all the possible sources in the whole Universe 
we should only replace dlogM/dE in the integral (jlip by 
the corresponding expression for photons, Eq. (fTUI) . and 
the cosmic-ray luminosity L by the photon luminosity 
L-v = ilL, Eq. (El). The order-of-magnitude estimate then 



, 7 - T)L, Eq 
reads 



10 D 



eV 



s • sr 



This total gamma-ray emission associated with particle 
acceleration in all the sources does not exceed the ob- 
served value of the diffuse gamma-ray background [4l| . 

We have also performed a more detailed numerical sim- 
ulation of the secondary gamma-ray flux. The injection 
spectrum of the curvature photons is similar to the syn- 
chrotron one WM , 
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FIG. 4. Gamma-ray fluxes predicted by the same models as 
shown in Fig. [2] 



18 1n 20 



POO 

w 2V3Z 2 e 2 j— / K 5/3 (x)dx, 



where 



Wc = 3- 

is the "critical" frequency (for higher frequencies, radi- 
ation is negligible), u is the frequency of the radiated 
photons, K 5 / 3 is the Macdonald function, 7 = (£/m) 
is the particle Lorentz factor, R is the curvature radius 
of the particle trajectory (for curvature radiation, R is 
constant and it is equal to the curvature radius of the 
magnetic lines). The total flux of the photons from a 
single source is related to the cosmic-ray flux by the co- 
efficient 77, Eq. (EJ). For the best-fit spectrum of Sec. MI Bl 
we used the same code to describe the propagation of the 
accompanying gamma rays and to calculate the observed 
gamma-ray flux. The result is presented in Fig. [4] One 
may see that the diffuse gamma-ray upper limit is satis- 
fied. 



V. CONCLUSIONS AND DISCUSSION 

We have constructed and studied a toy model of 
UHECR acceleration in the vicinity of numerous and var- 
ious supermassive black holes in centers of galaxies. The 
model assumes that: 
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• cosmic-ray particles are accelerated by the regular 
electric field within a few R$ from the SMBH 23] ; 
the field configuration is given by the solutions 
of Refs. [24], and is fully determined by the 
SMBH mass M , its angular momentum a and the 
magnetic- field normalization Bo; 

• all cosmic-ray particles accelerated near a given 
SMBH have similar initial conditions and there- 
fore all are accelerated up to one and the same 
energy limited by the curvature-radiation losses; 
this maximal energy is calculated in the model 
and depends on M only (provided Bq is a given, 
model-dependent, function of M ; dependence from 
a smooths this monochromatic spectrum insignifi- 
cantly) ; 

• the mean flux of a source (which accounts for the 
fraction of the sources where this mechanism does 
work) depends from M in a power-like manner; the 
normalization and the exponent are two free pa- 
rameters of the model (the best fit to the cosmic-ray 
spectrum indicates that this dependence is weak); 

• the concentration of sources is determined by 
the redshift-dependent SMBH mass function taken 
from astrophysical literature. 

Within these assumptions and given the Bq(M) relation 
is fixed (we considered three popular choices for it), the 
model has two free parameters which we find by fitting 
the cosmic-ray spectrum at the Earth to the experimental 
data. With parameters fixed in this way, we subject the 
model to several further tests which it passes succesfully: 

1 . the concentration of sources is large enough to sat- 
isfy the constraints from absence of clustering in 
UHECR arrival directions; 

2. the luminosity of a particular source, determined 
by the flux normalization and concentration, is not 
too high; 



3. secondary gamma rays from distant sources do 
not overshy the measured GeV diffuse gamma-ray 
background. 

Given the success of the toy model, it is interesting to dis- 
cuss its possible refinements. The assumptions we have 
made within the model are quite robust and realistic. 
One subtle point is related to the value of the SMBH 
angular momentum a which may vary from one black 
hole to another. However, these variations are probably 
modest given the scaling relation between the mass and 
the angular momentum of cosmic black holes proposed 
in Ref. (43|. The precision of predictions may be im- 
proved with more realistic modelling of the acceleration 
mechanism, in particular, with account of the charge con- 
centration at the SMBH and in the acceleration region, 
of the finite thickness of the accretion disk etc. Ulti- 
mately, this approach might give answer to the question, 
which particular SMBHs are strong sources and which 
are not, thus determining the (presently free) parameter 
theoretically. However, this is a complicated task and is 
far beyond the scope of the present work. 

One possible question concerns the low-energy part of 
the spectrum where, as is clearly seen from Fig. |U the 
contribution of the mechanism we discuss is insufficient 
to explain the observed spectrum due to the depletion of 
the SMBH mass function at low masses. It is tempting 
to speculate that this depletion is compensated by a huge 
contribution of the SMBH in our own Galaxy which, in- 
deed, has the appropriate mass. A quantitative analysis 
of this proposal requires, however, a much more precise 
study of physical properties and possibilities for particle 
acceleration close to the Galactic Center. 
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